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ABSTRACT 

H2 formation in metal-free gas occurs via the intermediate H" 



or H, ions. Destruction 



of these ions by photodissociation therefore serves to suppress H2 formation. In this 
paper, I highlight the fact that several processes that occur in ionized primordial gas 
produce photons energetic enough to photodissociate H~ or and outline how to 
compute the photodissociation rates produced by a particular distribution of ionized 
gas. I also show that there are circumstances of interest, such as during the growth of 
H 11 regions around the first stars, in which this previously overlooked form of radiative 
feedback is of considerable importance. 
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1 INTRODUCTION 

It has long been known that ro-vibrational line emission 
from molecular hydrogen, H2, is the dominant cooling pro- 
cess in gas of primordial composition at temperatures 200 < 
T < 10 4 K ilSaslaw fc Zipovl Il967l; iPeebles fc Dickd ll96S " 



iMatsuda. Sato, fc Takedall 19691 : Iffirasawa Ill969h . Moreover, 
recent work has made it clear that H2 cooling also plays 
a key role in the ev olution of low-density metal-poor gas 
ijjappsen et al.l 2007). In primordial gas and in low density 
metal-poor gas, H2 formation on dust is ineffecti ve, and H2 
form s primarily through the gas-phase reactions (|McDowelll 
ll96ll ; |Peebles fc Dickdll96sl ) 



H + e" 
H~ +H 



H-+7: 

H 2 + e~, 



(1) 

(2) 



with a smaller contribution coming from (|Saslaw fc Zipovl 

Gas3) 



H + H+ 
H++H 



HJ + 7, 
H2 + H 4 



(3) 
(4) 



Although only a small fractional abundance of H2, of order 



10" 



is produced by these reactions, many authors have 



shown that this is more than sufficient to allow the gas to 
cool within a Hubble time, and to allow star formation to oc- 
cur ( s ee e.g. iTegmark et afll 19971: iBromm. Coppi. fc Larsonl 



12002) : lAbel Bryan, fc Normanll2002l ; lYoshida et allbOOoC 

Once stars have formed, however, their radiation 
can interfere with the gas-phase production of H2. Most 
of the attention devoted to the study of so-called "ra- 
diative feedback" has focused on the ultraviolet radi- 
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ation from massive stars, which can suppress further 
star forma tion by photoionizing atomi c and molecular 



hydr o gen llHaiman. Thoul fc Loebl 19961: Kitavama et al 



2004 IWhalen. Abel fc Norman! 12004 IWhalen fc Norman 



20061 : lAbel. Wise fc Brvanl 120071) and by photodissociat- 
ing H? llHaiman. Rees fc Loeblll997l: lOmukai fc Nishi|[l999l : 
iHaiman. Abel fc Reesl l2000l : Idover fc Brandl l200ll ). with 
the latter process able to operate at far greater distances 
from the star, owing to the low opacity of the gas at fre- 
quencies redwards of the Lyman limit. On the other hand, 
radiative feedback due to the photodissociation of H~ and 
H+ 



H" +7 
Hj + 7 



H + e" 
H + H H 



(5) 
(6) 



has attracted little direct study. One reason for this is 
the speed of reactions [5] and [3] both have rate coefficients 
k ~ 10 -9 cm 3 s _1 , and so for reactions \5\ or [6] to be ef- 
fective at suppressing H2 formation, they must occur at a 
rapid rate. The other main reason for the comparative ne- 
glect of these processes is the nature of the radiation sources 
typically assumed in studies of radiative feedback. Work to 
date has focused on feedback from massive stars or from 
AGN. These sources are bright in the far-UV, and previous 
studies have found that feedback due to H2 photodissoci- 
ation by UV photons from these sources becomes effective 
long before feedback due to H~ or Hj photodissociation by 
optical or near-IR photons from the same sources (see e.g. 
iMachacek. Bryan fc AbelbOOlT ). 

However, one issue that does not appear to have been 
previously considered in any great detail is that fact that 
once sources of ionization such as massive stars or AGN ex- 
ist, the ionized gas that they produce will act as a secondary 
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source of radiation. Radiative feedback due to the pho- 
todissociation of H2 by Lyman- Werner ban d photons pro- 
duced b y recombining He + was considered bv lYoshida et al] 
(2006b) and found to be unimportant. However, radiative 
feedback due to the photodissociation of H" or Hj by re- 
combination emission or free-free emission from the ionized 
gas does not appear to have been considered prior to now. 

In this paper, I investigate this form of feedback and 
show that there are circumstances of cosmological impor- 
tance in which it can be highly effective. The calculations 
presented in this paper involve a number of simplifications, 
so as to facilitate the exploration of a variety of different 
scenarios with minimal computational effort. Nevertheless, 
they should produce results that are correct to within an or- 
der of magnitude, and so serve to highlight which scenarios 
require more detailed study. 

The layout of this paper is as follows. In f]2] I outline the 
method used to compute the photodissociation rates pro- 
duced by a given distribution of ionized gas. In Sj3] I use 
this method to investigate the effects of H~ and H^~ pho- 
todissociation in the neutral gas surrounding an expanding 
H n region created by a population III star ( i|3.ip . in a neu- 
tral cloud embedded in a large H 11 region f £|3. 2\ and in the 
centre of a recombining 'fossil' H 11 region ( N3.3|l . Finally, in 
ij4] I summarize what these results tell us about the scenarios 
in which H~ and Hj photodissociation is important. 



2 MODELLING EMISSION FROM IONIZED 
GAS 

The radiation flux at a location x produced by emission from 
an ionized volume V is given by: 



F(v) 



1 

4% 



-r(„,x',x) n c( x ') 7i(^; x>i(x') 

|x' -x| 2 



dx', (7) 



where t(u, x', x) is the optical depth at frequency v between 
x and x', n c is the number density of electrons, m is the num- 
ber density of ions of species i, ji is the emission coefficient 
for species i (with units erg cm 3 s _1 Hz -1 ), and where we 
sum over all ionic species present in the gas. For clarity, in 
the simplified examples presented in this paper, I restrict my 
attention to emission from recombining H + ions. In normal 
circumstances this will produce the bulk of the emission, 
with a small additional contribution coming from helium, 
and with only negligible amounts of radiation coming from 
other ions. This simplification allows us to rewrite equation]?] 
as: 



-t(«,x / ,x) n c (x')n H+ (x ; )7 H+ (v; x') 



dx', (8) 



where n H + is the number density of H + ions and where 
7 H + is the emission coefficient for emission from H + . We 
can further simplify this equation by making the assumption 
that the temperature of the ionized gas is uniform. In this 
case 7 H + (u) becomes independent of position and can be 
moved outside the integral: 



F(u) 



7h+ W f „-t(„,x',x) w c (x')n H +(x') 



47T 



dx'. 



(9) 



A final simplification comes if we assume that the gas is 
optically thin. At the densities of interest in this work, 



this is a good approximation for most of the emission, 
as the continuum opacity of metal- free gas is very small 
l|Lenzuni. Chernoff fc Salpeterl Il99lf ). The exceptions are 
the bound-free emission produced by recombination directly 
into the n = 1 ground state, and the Lyman series lines. As 
I discuss at greater length in the next section, this emission 
can be taken into account by assuming that case B recom- 
bination applies. 

The assumption of optically thin gas allows us to rewrite 
Equation [9] as: 



where 



4tt 



n e (x')n H +(x') 



Ix' 



dx'. 



(10) 



(11) 



With these simplifications we have reduced the problem of 
computing the flux into two separate, simpler problems: that 
of computing 7 H + (y) and that of evaluating the integral Iv- 
For a given distribution of ionized gas, numerical evaluation 
of the integral is trivial, and so the only real remaining dif- 
ficulty is the evaluation of 7 H + (u). I discuss this in the next 
section. 



2.1 Computing the emission coefficient 

2.1.1 Bound-bound emission 

In my treatment of the effects of recombination, I assume 
that Case B applies. In the classical Case B, the gas is op- 
tically thick in both the Lyman continuum and the Lyman 
series lines. Therefore, emission produced by recombination 
directly to the n = 1 ground state will be immediately re- 
absorbed, with the result that only recombinations to ex- 
cited states with n ^ 2 actually result in a net decrease in 
the number of ionized hydrogen atoms. Every recombination 
to a state with n ^ 2 will be followed by a radiative cascade 
as the atom attempts to reach the ground state and some 
fraction of the photons produced during these cascades will 
be capable of photodissociating H~ or Hj. 

The photodissociation threshold of H~ is 0.755 eV, and 
so H~ can be photodissociated by any photons in the Lyman 
series or Balmer series, and by most photons in the Paschen 
series (starting with Paschen-/?). H^", with its larger pho- 
todissociation threshold of 2.65 eV, can be dissociated by 
any Lyman series photon, and by most Balmer series pho- 
tons (starting with H7), but is not affected by Paschen series 
photons. 

For the emissivities of the lines in the Balmer and 
Paschen series, I use the values computed for Case B by 
IStorev fc Hummer! l|l995h , who tabulate frequency-averaged 
emissivities for transitions from all excited states up to 
n — 50. As the widths of these emission lines are much 
smaller than the frequency range over which the photodis- 
sociation cross-sections of H~ or Hj vary substantially, ac- 
curate modelling of the line profiles is unnecessary and so 
for simplicity I model the lines as delta functions. 

In Case B, all Lyman series photons more energetic 
than Lyman-a are eventually degraded to Lyman-a pho- 
tons pl us one or more lower ene rgy photons. While the re- 
sults of lStorev fc Hummer! l|l995h include these lower energy 
photons, they do not account for the eventual fate of the 
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Lyman-a photons. In the classical Case B, these photons 
can never escape from the ionized region: the optical depth 
of the Lyman series lines is assumed to be infinite. In reality, 
the optical depth is finite and some photons can escape from 
the gas, although generally not before they have scattered 
multiple times. There is also a small chance that an excited 
hydrogen atom in the 2p state will reach the ground state 
not by direct radiative decay to the Is state, but rather by 
undergoing a collisional transition to the 2s state, followed 
by two-photon decay to the ground state. Over time, this 
process results in the net loss of a large number of Lyman-a 
photons from the ionized gas. 

The relative importance of these two processes de- 
pends upon a number of factors, such as the density of the 
H n region and the Lyman-a optical depth of the gas. For 
simplicity, I consider in this paper only the limiting case in 
which two-photon decay dominates. In this case, the contri- 
bution to 7 H + (v) arising from two-phot on emission can be 
written as dFernandez fc Komatsu|[200r3) 



72phM = -^-P{v/v Lya )a B , 



(12) 



where v^ ya is the frequency at line center of the Lyman- 
a line, ob is the Case B recombination coefficient and 
P(y)dy is the normalized probability per two-photon de- 
cay of getting one photon in the inter val Ay = Av/vhyg. 
For q b, I use the values tabulated by IStorev fc Hummer] 
(Il995h. while for P(y) I use a fitting function taken from 
iFernandez fc Komatsul i|2006l ): 



P(y) = 1.307- 



where 



2. 6272 2 + 2.563/ 



0.5. 



51.692 



(13) 



2.1.2 Bound-free and free- free emission 

The contribution to 7 H+ (i/ ) from bound- free and free- f ree 
emission can be written as (|Fernandez fc Komatsul [20061 ) 



7bf+ff 



6.84 x 10 _3i 
T 1 / 2 



-his/kT 



_ \~y x n e " 
gs + > J flfb(n) 



n=2 



(14) 



where T is the gas temperature, x n — lRyd/(fcTn 2 ), gg is 
the thermally averaged Gaunt factor for free-free emission 
and gib{n) is the Gaunt factor for free-bound emission from 
recombination into level n. The value of this expression has 
be en tabulated as a function of frequency and temperature 
bv lFerlandl l|l980h . In Figure [T] I plot the value of 7bf+ff as 
a function of photon energy for photons with 0.755 < hv < 
13.6eV, assuming a gas temperature T = 10 4 K. In the same 
figure, I also plot the two-photon contribution 72 P h- At low 
energies, close to the H~ photodissociation threshold, the 
bound-free and free-free contributions dominate, while at 
high energies, the two-photon contribution dominates. It is 
also plain from the figure that if we were to disregard the 
two-photon emission, and to instead assume that all of the 
Lyman-a photons produced by the ionized gas eventually 
escape by scattering into the wings of the line, we would 
nevertheless obtain photodissociation rates for H~ and 
of the same order of magnitude as those derived below. 
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where the photodissociation cross section is (|Wisharj|l979l ) 



<r H -(y) = 7.928 x 10 J 



R 



ffH+(f)F(l/) 



I IV 



-d/A 



cr H + (^)7h+ {v) 

Iv I — - — r d^i 

hv 
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Figure 1. Emission coefficients, plotted as a function of pho- 
ton energy, for the bound-free and free-free emission (summed 
and plotted as solid line) and two-photon emission (dashed line) 
produced by an ionized gas of pure hydrogen with temperature 
T = 10 4 K. 



2.2 Evaluating the photodissociation rates 

Given 7 H + , it is easy to compute the photodissociation rates 
of H~ and H^" as functions of the volume integral Iv- The 
H~ photodissociation rate can be written as 



(15) 



(16) 



(17) 



at frequencies above the H~ photodissociation threshold 
hvth = 0.755 eV. The photodissociation rate is given 
by a similar expression 



(18) 



(19) 



where the photodissociation cross section at photon energies 
11.27 eV > hv > hv tb is 

ct h+ (v) = dex [-40.97 + 82.01?? - 93.22?7 2 + 34.89r? 3 ] , (20) 

and at photon energies hv > 11.27eV is 

a H + {v) = dex [-30.26 + 37.947? - 34.03?? 2 + 8.89T/ 3 ] , (21) 

where r/ is the photon energy in Rydbergs, and where 
the photo dissociation threshold energy hv tb — 2.65 eV = 
0.195 Ryd (|Dunnlll968l ). 
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H II region temperature (K) 







5000 


10000 


20000 




Continuum 


3.49 


1.78 


1.01 




Line 


23.8 


3.44 


0.720 




Total 


27.3 


5.21 


1.73 




Continuum 


0.100 


0.048 


0.024 




Line 


1.0 x 10- 10 


3.0 x 10- 11 


9.7 x 10~ 12 




Total 


0.100 


0.048 


0.024 


Table 1. Value of R pd 


i/Ivi m units of 10 29 cm 5 s 


\ fori = H" 



tures. The size of the contributions from continuum emission and 
line emission is noted. 



The values of the frequency integrals in equations [TS] 
and [19] are independent of the spatial distribution of the 
ionized gas, but do depend on its temperature. For gas tem- 
peratures in the range 3000 < T < 30000 K, their values are 
fit to within 5% by the expressions: 



R 



pd.H- 



Iv 
and 
R 



dex [- 



28.28 - 2.04 U + 1-28*4 



= dex [-30.33- 1.01*4 + 0.24*1] cm 5 s~ 



(22) 



(23) 



where *4 = log(T/10 4 K). Table [T] shows the relative size 
of the line and continuum contributions for three different 
gas temperatures: T = 5000 K, 10 4 K and 2 x 10 4 K. Two 
main points are worthy of note. Firstly, there is an obvi- 
ous temperature dependence: the hotter the Hn region, the 
lower the photodissociative flux. Secondly, line emission is 
far more important in the case of H~ photodissociation that 
in the case of Hj photodissociation. In the former case, line 
emission is the dominant contribution when the tempera- 
ture of the ionized gas is small and is significant even when 
T is large. In the latter case, the contribution from line emis- 
sion is always negligible. This difference is a consequence of 
the fact that the Hj photodissociation cross-section is very 
small compared to the H~ photodissociation cross-section 
at the energies which correspond to Balmer series photons. 
The main contribution to photodissociation comes from 
photons with energies ~ 6 eV or above and hence in this 
case the continuum processes dominate. 



3 RESULTS 

In order to establish whether or not the photodissociaton of 
H~ and by emission from ionized hydrogen is ever an 
important process, I compute in this section the amount of 
flux produced in several scenarios of cosmological interest, 
along with the size of the resulting photodissociation rates. 
As my aim in this paper is simply to produce estimates of 
these values with the correct order of magnitude, the model 
H II regions considered here are all highly simplified com- 
pared to realistic systems. Nevertheless, the results of these 
simple models should be a good guide as to whether more 
detailed numerical modelling is justified. 



3.1 Expanding H II region 

The first simple scenario considered here involves emission 
from an expanding Hn region within a small protogalac- 
tic halo. Because the first stars to form in the Universe 
are predicted to be very massive llAbel, Bryan, fc Normanl 
12002 " 



p redicted to be very mass ive (Abel, B ryan. & IN orman 
alBromm. Coppi. fc LarsorJl2002t lYoshida et aiTl2006a: 



iGao et alJ 2006) and theref ore to be emitte rs of large num- 
bers of ionizing photons |Schaererl |2002h . this scenario 
has been studied numerically by a number of authors 
llWhalen. Abel fc Normanl 120041: [Whalen fc Normanl 120061: 
Alvarez. Bromm fc Shapiroll2006|: Johnson. Greif fc Bromml 



20061 : ISusa fc Umemural 120061 : lAbel. Wise fc Brvanl 12007m 



This work has given us a good understanding of the basic 
chain of events. The nascent H 11 region initially grows very 
rapidly, but recombinations in the dense gas surrounding the 
ionizing source quickly cause its growth to slow. After less 
than 10 s yr (assuming a 100 Mq star), the ionization front 
bounding the H 11 region has become a slow, D-type front 
that drives a shock ahead of itself into the dense neutral 
gas. At radii greater than the shock radius rs, the density 
profile of the gas remains almost the same as before the 
switch-on of the ionizing source. Within the H 11 region, on 
the other hand, the density profile is flat, indicating that 
the density of the ionized gas is almost constant. At this 
point, the evolution of the H 11 region is dominated by the 
pressure-driven expansion of the ionized gas: as the gas ex- 
pands, the recombination rate falls, and so more ionizing 
photons become available for ionizing previous neutral gas, 
expanding the amount of mass contained in the H 11 region. 
If the density profile of the protogalactic gas is steep enough 
- it must fall off with radius more quickly th an p <x r _3//2 
l|Franco. Tenorio-Tagle. fc Bodenheimerlll990l ) - and if the 
massive star shines for long enough, then eventually the ion- 
ization front will undergo a transition back to a supersonic 
R-type front, and will soon thereafter break out of the proto- 
galaxy. However, even if the source switches off before break- 
out occurs, most of the ionized gas will nevertheless escape 
from the protogalaxy, as the mean outward radial velocity 
of this material is typically several times higher than the 
escape velocity of the halo. 

Although the strong Lyman-Werner band emission 
from the massive star will quickl y dissociate diffuse H2 
throughout much of the halo (|Omukai fc NishU 1 19991 : 
iGlover fc Brand! l200ll ). the presence of abundant free elec- 
trons in the partially ionized region ahead of the ioniza- 
tion front can ca t alyze H2 formation in this region, and 
ISusa fc Umemural (2006) have shown that the shielding pro- 
vided by this H2 can allow H2 to survive in dense cores else- 
where in the protogalaxy. However, these calculations did 
not take into account the possible effects of recombination 
emission from the ionized gas. 

To compute the emission from a typical protogalactic 
H 11 region, we assume, for simplicity, that it is spherically 
symmetric, with radius Ri, and that the gas within it has 
a uniform density ni. Inspection of the results from de- 
tailed three-dimensional numerical simulations of H 11 region 
growth in these protogalactic systems suggests that these 
should be reasonable approximations. In that case, the value 
of Iv at a distance D > Ri from the center of the H 11 region 
is given by 



Hj Rif(xv) 



(24) 
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Figure 2. Value of the function /(id), defined in Equation 1251 
for a wide range of xv 



where 



/(zd) = 



1 + 



(iD - X) 2 



Ax, 



(25) 



and where id = D/R\. In deriving this expression I have 
assumed that the gas is fully ionized, with n H + — n c = m. 
The value of f(xu) is plotted in Figure [2] as a function of 
id- In the limit that Xu ~* 0, we see that Jy — > n 2 Rj, 
while for id = 1, iV = 0.5n 2 iii and in the limit of large D, 
Iv ^ (l/3D 2 )nfRi. 

If we assume that the temperature of the ionized gas is 
2 x 10 4 K, then at a point just outside the Hn region, the 
H~ photodissociation rate is given by 



-ii / m \ 2 ( _Rr_\ _i 
Vl cm" 3 / \ 1 pc ) S 



-Rpd.H- = 2.67 x 10 



and the Hj photodissociation rate is given by 



R . „+ = 3.76 x 10" 

pd,H, 



13 ( ni \ 2 ( R x \ _i 
Ucm- 3 / \ lpc ) S 



(26) 



(27) 



For comparison, the rate at which H is destroyed by asso- 
ciative detachment (reaction [3J is given by 



R cH - = 1.3 x 10~ 9 n H s~\ 



(28) 



where nn is the number density of atomic hydrogen, and 
where we have a dopted the rate coefficient determined for 
this reaction by ISchmeltekopf et al.l |l967). As noted by 
iGlover. Savin fc Jappsenl (|2006l ). the rate of this reaction 
is uncertain by almost an order of magnitude, but for the 
initial order-of-magnitude study given in this paper, the 
use of a value from close to the middle of the range of 
uncertainty is probably justified. The rate at which 
is destroyed by charge transfer ( reaction [4) is given by 
jKarpas. Anicich fc Huntressll 19791 ) 



6.4 x 10~ 10 tc h s~\ 



(29) 



In this particular example, the density of atomic gas 
just outside the H II region will be similar to the density of 
the ionized gas within the H n region, implying that nn ~ m. 



Therefore, H will be destroyed primarily by photodissoci- 
ation, rather than associative detachment, if 



ni >50 



Rx 
1 pc 



(30) 



Similarly, photodissociation will dominate the destruction of 



HJ if 



Ri 



m > 2400 [ | cu) 

lpc 



(31) 



How do these numbers compare with the results of the nu- 
merical simulations of H n region growth mentioned above? 
To take o ne particular example , consi der the H n region sim- 
ulated bv lAbel. Wise fc Brvar] \200H ). At t = 0.1 Myr after 
the switch-on of the central ionizing source, the H II region in 
their simulation has a mean density n\ ~ 500cm~ 3 and a size 
J2i ~ 1 pc. At a later time, t = 1 Myr, its size has grown to 
Ri ~ 20pc, while its mean density has fallen to n\ ~ 10cm -3 . 
In both cases, the diffuse flux from the Hn region is strong 
enough to dominate the destruction of H~ in the surround- 
ing gas. It does not dominate the destruction of H^", but is 
important at the 10-20% level. Only once t = 2.7 Myr (at 
which point Ri ~ 50 pc and m ~ 1 cm -3 ) does the diffuse 
flux decrease to a level at which it no longer dominates the 
destruction of H~. However, by this point in the evolution of 
the protogalaxy, the expansion of the central H II region has 
accelerated almost all of the surrounding gas to velocities 
greater than the protogalactic escape velocity. Therefore, 
any H2 that does manage to form via H~ at t > 2.7Myr will 
not be retained by the protogalaxy. 

Although the preceding analysis is for gas located just 
outside the H 11 region, H~ photodissociation will actually 
dominate out to much larger distances. This can be easily 
seen if we compare the dependence on radius of -R p d,H- an d 
_R C H - . In the limit of large distances, -R p d,H- scales with the 
distance D as -R p a,H- x -D -2 , while at smaller distances the 
scaling is less steep. If the density distribution of the neutral 
gas is given by n oc D~ a , then R c n- oc D~ a . Therefore, if 
a > 2, then R c K - will fall off more steeply with increasing 
radius than i?p dH -. As simulations of popu lation III star 
format ion generally find that a ~ 2.2 (see e.g. lYoshida et ail 
l2006al ). this implies that if H~ photodissociation dominates 
near the H 11 region, it will also dominate throughout the 
remainder of the protogalaxy. 



3.2 Neutral cloud embedded in large H 11 region 

The second scenario considered in this paper is that of a 
cloud of neutral hydrogen embedded in a large cosmologi- 
cal H 11 region. This arrangement of ionized and neutral gas 
may be encountered within the earliest protogalaxies after 
the formation of the first stars if any clumps of gas exist 
within these protogalaxies that are dense enough to survive 
the passage of the expanding ionization fr ont (I-front) with- 
out being imme diately photoevaporated (|Susa fc Umemural 
l2006l ; ISusall2007i ). On larger scales, this arrangement can also 
be used to represent primordial minihalos that have been en- 
veloped by a large intergalactic H 11 region, but that have not 
yet been completely photoionized or photoevaporate d (see 
e.g. Illiev. Shapiro fc Ragall2005l ; lAhn fc Shap"Irol 120071 1. 

For simplicity, I assume that both the neutral cloud and 
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the surrounding H II region are spherically symmetric, with 
radii R c i and Ri respectively. The cloud is assumed to have 
a power-law density profile 



(32) 



with a < 3, where n c i is the density of the cloud at R = R c i, 
which can be written in terms of the cloud mass M c \ as 



n c i 



Mci 3 - a 
m H 47rR 3 ' 



(33) 



where run is the mass of a hydrogen atom. At early times, 
the outer edge of the cloud will be photoionized, but will 
not yet have had time to dynamically respond to the in- 
creased pressure, and so will maintain the same density pro- 
file as the neutral gas. The density of ionized gas at radii 
i?ci > R > -Rci. i, where R c i,i is the innermost extent of 
the ionized region, is therefore given by Equation 1321 above, 
while at radii R R c i I take the ionized gas density, m, to 
be constant. Finally, it is necessary to specify the distance D 
of the centre of the cloud from the centre of the H II region. 

These assumptions allow us to write the volume integral 
Iv as 



2 d^cl 



2a 



Rci 
Rci,: 



+ n 2 Rif(x D )-n 2 R cU 



(34) 



where id = D/Ri and where /(id) is defined by Equa- 
tion [25] The first term in this expression corresponds to the 
flux produced by the ionized outer edge of the cloud, while 
the second and third terms, taken together, correspond to 
the flux from the surrounding constant density Hn region. 

We can explore the behaviour of this expression by con- 
sidering a concrete example. Let us suppose that D — 0.5-Ri, 
so that id = 0.5, and that -R c i,i = 0.5-R c i and a = 2. Then 
Iv becomes 



^riciRd + 0.91nj-Ri 



2 r> 

nj Rci. 



(35) 



If -Ri ^> -Rd, then the second term dominates, unless the 
density contrast between the cloud and the surrounding 
H II region is large. The first term dominates only if n c i > 
0.6(-Ri/-Rd) 1//2 wi, while the third term is always negligible 
in these conditions. The required density contrast can be 
significantly reduced if we steepen the density profile of the 
cloud by increasing a, or increase the fraction of it which is 
already ionized by decreasing R c i,i- Altering D, on the other 
hand, has little effect on the required contrast. 

To proceed with our example, let us suppose that the 
H II region has a size Ri = 3 kpc and numb er density ri\ — 
l(T 3 cm" 3 j Johnson. Greif fc Brorrrm] |2006). that the cloud 
has a virial radius -R c i = 100 pc, and that n c i = 0.2 cm" 3 . 
In that case, the flux from the ionized edge of the cloud 
dominates, Iv — 2.9 x 10 19 cm -5 , and the resulting H~ and 
rlj photodissociation rates are R p d,u- 
and R , „+ 

pd,Hj 



6.9 x 10 12 s 1 respectively. In comparison, 



R. 



1.04 x 10" 



1 and R 



at 



. c ,h — " — " — " C>H;| 

Rd, i. Therefore, in this particular example the flux is of lim- 
ited importance. If we were to assume that ,R c i,i = 0.05-R c i, 
rather than -R c i,i = 0.5-R c i, however, then the photodissoci- 
ation rates would increase by a factor of a thousand, while 



the collisional rates would increase by only a factor of a 
hundred. In that case, more H~ would be destroyed by pho- 
todissociation than by associative detachment. On the other 
hand, the photodissociation of would still be of limited 
importance. 

This simple example demonstrates that in this situa- 
tion, radiative feedback due to the emission from the ionized 
gas is most effective during the final stages of the photoion- 
ization of the neutral cloud, when the highest density ma- 
terial is being photoionized. As we expect dynamical effects 
to have beco me important by this s tage in the photoioniza- 
tion process (|Ahn fc Shapiro! |2007| ). it is difficult to assess 
the ultimate importance of the feedback without performing 
more detailed calculations. 



3.3 Recombining fossil H II region 

The final s cenario examined her e is that of a so-called 'fossil' 
H II region (|Oh fc Ha iman 2003]), i.e. an Hn region in which 
the ionized source has switched off. For simplicity, I consider 
the evolution of the gas at the center of a spherically sym- 
metric H II region with constant density n\ and radius R\. If 
the gas is initially fully ionized, then Iv at the moment of 
switch- off is given by 



Iv = n\R\. 



(36) 



If we again assume that the temperature of the ionized gas 
is 2 x 10 4 K, then the flux produced by the ionized gas gives 
rise to photodissociation rates for H~ and that are given 
by 



R r , 



R 



5.3 x 10" 



7.4 x 10" 



Vl cm" 3 / 

o / ni 
V 1 cm" 



Ri 
1 kpc 

Ri 
1 kpc 



(37) 



(38) 



As time passes, however, the gas in the Hn region will re- 
combine. Since Iv oc x 2 , where x is the fractional ionization 
of the gas, this means that both R p< j,h- an d R pd H + will 
both decrease with time. If the temperature of the ionized 
gas were to be kept fixed, then these rates would also fall off 
as x 2 . In fact, it is likely that some cooling of the gas will 
occur and so the effects of the decrease in x will be offset 
to some extent by an increase in the emission coefficient. 
Nevertheless, it is clear from Table [T] that changes in the 
emission coefficient will alter the rates by at most an order 
of magnitude, while a large decrease in x may alter them 
by many orders of magnitude. Therefore, consideration of 
the simplified case in which T is kept fixed serves to illus- 
trate the basic behaviour of the gas without requiring one 
to model its thermal evolution. 

For a intergalactic Hn region with the same size and 
density as in the previous section (i.e. Ri = 3 kpc and ni = 
10~ 3 cm -3 ), we therefore have 



R 
R 



pd,H- 



pd,H+ 



= 1.59 x 10 -13 i 2 s~\ 
= 2.22 x lO" 15 ! 2 s" 1 . 



(39) 
(40) 



The destruction rate of H by associative detachment in the 
same conditions is given approximately by 



Rr 



1.3 x 10" 12 (1 - x) s" 1 . 



(41) 
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Clearly, -R p d,H- > Rc,n- om y if (1 — x ) ^ !• However, if 
the fractional ionization is as high as this, then H~ will also 
be destroyed rapidly by mutual neutralization with protons 

rT+H+^H + H. (42) 

The rate coefficient of this reactio n is uncertain, possibly by 
as mu ch as an order of magnitude (jGlover. Savin fc Jappsenl 
2006). If we take the s mallest of the values quo ted in the 
literature, the rate of iDalgarno fc Leppl ([l987), then we 
find that the mutual neutralization rate in our example 
H n region is 

-Rmn.H- = 7X W~ W T- 1/2 X s"\ (43) 

This is comfortably larger than -R p( j,H- f° r an temperatures 
in our range of interest. Therefore, in this particular sce- 
nario, H~ photodissociation is unimportant. Comparison of 
the photodissociation rate of with the rate at which it 
is destroyed by charge transfer with H (reaction [4| or by 
dissociative recombination 

Ut + e~ -» H + H, (44) 

leads to a similar conclusion regarding Hj. 

Can we avoid these conclusions by considering a larger 
H II region? If we suppose that we have somehow managed 
to create a fossil Hll region with Ri = 3 Mpc, rather than 

3 kpc, then it is easy to see that the resulting photodisso- 
ciation rates would be a thousand times larger than in the 
case considered above. Even so, associative detachment still 
dominates the destruction of H~ if x < 0.1. Similarly, pho- 
todissociation is unimportant in comparison to dissociative 
recombination or charge transfer for x < 0.5. Since H2 for- 
mation in cooling, rec ombining gas become s significant only 
once x < 0.1 (see e.g. lOh fc Haimanll2002h . it is clear that 
even in this somewhat unrealistic case, photodissociation of 
H~ and is unimportant. 

The other way in which to increase the importance of 
photodissociation is by increasing the density of the ionized 
gas. If, instead of a large intergalactic H 11 region, we consider 
a small interstellar H 11 region with m = 100 cm -3 and Ri = 
10 pc, then we find that 

Rpd,H- = 5.3 x 10 - Vs -1 , (45) 
^pd.Hj = 7.4xlO _ Vs _1 . (46) 

In this case, R p dji- > ^oH- omv if ar > 0.15, while 
R pd H + > R c H + only if x > 0.59, so again the effects of 
photodissociation are unimportant. 

4 CONCLUSIONS 

The results from the simple models considered in the previ- 
ous section demonstrate that in some circumstances, bound- 
free, free-free and two-photon emission from ionized gas can 
combine to produce a significant H~ photodissociation rate. 
The effect is very sensitive to the density of the ionized gas, 
since both the recombination rate and the free-free emission 
rate scale as the square of the density. However, the required 
densities are not excessively high and there are a number of 
situations in which we may encounter them. 

The particular examples considered here in which H~ 
photodissociation proves to be important are the early 



growth of an H 11 region around a population III star, prior to 
its break-out from the confining protogalaxy f t]3.1[) and the 
late stages of the photoionization of a neutral cloud (or mini- 
halo) by an external radiation source f £|3.2[) . On the other 
hand, photodissociation of H~ is not important within so- 
called 'fossil' H 11 regions, owing to the rapid decrease of the 
photodissociating flux associated with the decrease in the 
fractional ionization of the recombining gas. 

The destruction of H \ by emission from the ionized gas 
proves to be a less important effect, due to the much smaller 
photodissociation rate (compared to H~) that results from 
the same amount of emission. The main reason for this is 
that the photodissociaton cross-section is small (or zero) 
at the photon energies where much of the energy from the 
ionized gas is radiated, becoming significant only for pho- 
tons with energies greater than about 6 eV. Nevertheless, 
even if H^~ survives where H~ does not, the net effect is still 
a significant reduction in the H2 formation rate, as H2 for- 
mation via Hj~ occur s much more slowly than H2 formation 
via H" (|Gloverll2003r i. 



ACKNOWLEDGMENTS 

The author would like to thank T. Greif and the anonymous 
referee for their comments on an earlier draft of this paper. 



REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 
295, 93 

Abel, T., Wise, J. H., Bryan, G. L. 2007, ApJL, in press; 
|astro-ph/0606019| 

Ahn, K., & Shapiro, P. R. 2007, MNRAS, 375, 881 
Alvarez, M. A., Bromm, V., & Shapiro, P. R. 2006, ApJ, 
639, 621 

Bromm, V., Coppi, P. S., & Larson, R. B. 2002, ApJ, 564, 
23 

Dalgarno, A., & Lepp, S. 1987, in Astrochemistry, ed. M. S. 

Vardya & S. P. Tarafdar (Dordrecht: Reidel), 109 
Dunn, G. H. 1968, Phys. Rev., 172, 1 
Ferland, G. J. 1980, PASP, 92, 596 
Fernandez, E. R., & Komatsu, E. 2006, ApJ, 646, 703 
Franco, J., Tenorio-Tagle, G., & Bodenheimer, P. 1990, 

ApJ, 349, 126 

Gao, L., Abel, T., Frenk, C. S., Jenkins, A., Springel, V., 

& Yoshida, N. 2006, |astro-ph/0610174| 
Glover, S. C. O., & Brand, P. W. J. L. 2001, MNRAS, 321, 

385 

Glover, S. C. O. 2003, ApJ, 584, 331 

Glover, S. C. O., Savin, D. W., & Jappsen, A.-K. 2006, 
ApJ, 640, 553 

Haiman, Z., Abel, T., & Rees, M. J. 2000, ApJ, 534, 11 
Haiman, Z., Rees, M. J.,& Loeb, A. 1997, ApJ, 476, 458 

(Errata in ApJ, 484, 985) 
Haiman, Z., Thoul, A., & Loeb, A. 1996, ApJ, 464, 523. 
Hirasawa, T., 1969, Prog. Theor. Phys., 42, 523 
Hummer, D. G., & Storey, P. J. 1987, MNRAS, 224, 801 
Iliev, I., Shapiro, P. R., & Raga, A. C. 2005, MNRAS, 361, 

405 



8 S. C. O. Glover 



Jappsen, A.-K., Glover, S. C. O., Klessen, R. S., & Mac 
Low, M.-M. 2007, ApJ, in press; |astro-ph/0511400| 

Johnson, J. L., G reif, T. H., & Bromm, V. 2006, 
|astro-ph/0612254| 

Karpas, Z., Anicich, V., & Huntress, W. T. 1979, J. Chem. 
Phys, 70, 2877 

Kitayama, T., Yoshida, N., Susa, H., & Umemura, M. 2004, 
ApJ, 613, 631 

Lenzuni, P., Chernoff, D. F., & Salpeter, E. E. 1991, ApJS, 
76, 759 

Machacek, M. E., Bryan, G. L., & Abel, T. 2001, ApJ, 548, 
509 

Matsuda, T., Sato, H., & Takeda, H. 1969, Prog. Theor. 

Phys., 42, 219 
McDowell, M. R. C. 1961, Observatory, 81, 240 
Oh, S. P., & Haiman, Z. 2002, ApJ, 569, 558 
Oh, S. P., & Haiman, Z. 2003, MNRAS, 346, 456 
Omukai, K., & Nishi, R. ApJ, 518, 64. 
Peebles, P. J. E., & Dicke, R. H. 1968, ApJ, 154, 891 
Saslaw, W. C, & Zipoy, D. 1967, Science, 216, 976 
Schaerer, D. 2002, A&A, 382, 28 

Schmeltekopf, A. L., Fehsenfeld, F. C, & Ferguson, E. E. 

1967, ApJ, 118, L155 
Storey, P. J., Hummer, D. G. 1988, MNRAS, 231, 1139 
Storey, P. J., Hummer, D. G. 1995, MNRAS, 272, 41 
Susa, H., & U memura, M. 2006 ApJ, 645, L93 
Susa, H. 2007, |astro-ph/0701172| 

Tegmark, M., Silk, J., Rees, M., Blanchard, A., Abel, T., 

& Palla, F. 1997, ApJ, 474, 1 
Whalen, D., Abel, T., & Norman, M. L. 2004, ApJ, 610, 

14 

Whalen, D., Norman, M. L., 2006, ApJS, 162, 281 
Wishart, A. W. 1979, MNRAS, 187, 59P 
Yoshida, N., Abel, T., Hernquist, L., & Sugiyama, N. 2003, 
ApJ, 592, 645 

Yoshida, N., Omuaki, K., Hernquist, L. & Abel, T. 2006a, 
ApJ, 652, 6 

Yoshida, N., Oh, S. P., Kitayama, T., & Hernquist, L. 
2006b, ApJ, in press; |astro-ph/0610819| 



